function domfront(p, e)
% domfront(p, e)

p1 = [[p(1); p(1); p(1); p(1)], [p(2); p(2); p(2)+e(2); p(2)+e(2)], [p(3); p(3)+e(3); p(3)+e(3); p(3)], [0.25;0.25;0.25;0.25]];
p2 = [[p(1); p(1); p(1)+e(1); p(1)+e(1)], [p(2); p(2); p(2); p(2)], [p(3); p(3)+e(3); p(3)+e(3); p(3)], [0.5;0.5;0.5;0.5]];
p3 = [[p(1); p(1)+e(1); p(1)+e(1); p(1)], [p(2); p(2); p(2)+e(2); p(2)+e(2)], [p(3); p(3); p(3); p(3)], [0.75; 0.75; 0.75; 0.75]];

%caxis([0 1])

m = colormap;

%h = patch([p1(:, 1), p2(:, 1), p3(:, 1)], [p1(:, 2), p2(:, 2), p3(:, 2)], [p1(:, 3), p2(:, 3), p3(:, 3)], [p1(:, 4), p2(:, 4), p3(:, 4)]);
h = patch([p1(:, 1)], [p1(:, 2)], [p1(:, 3)], [p1(:, 4)]);
set(h, 'EdgeColor', m(17, :));
%set(h, 'EdgeColor', 'none');
h = patch([p2(:, 1)], [p2(:, 2)], [p2(:, 3)], [p2(:, 4)]);
set(h, 'EdgeColor', m(33, :));
%set(h, 'EdgeColor', 'none');
h = patch([p3(:, 1)], [p3(:, 2)], [p3(:, 3)], [p3(:, 4)]);
set(h, 'EdgeColor', 'none');
%get(h)
